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Abstract 

Resonant formation of the muonic molecule dt^i in t^i atom collision with condensed H/D/T 
targets is considered. A specific resonance correlation function, which is a generalization of the 
Van Hove single-particle correlation function, is introduced to calculate the resonant-formation 
rate in such targets. This function is derived in the case of a polycrystalline harmonic solid. 
Also is found a general asymptotic form of the resonance correlation function for high momentum 
transfers, valid for any solid or dense-fluid hydrogen-isotope target. 

Numerical calculations of the rates are performed for solid molecular hydrogens at zero pressure, 
using the Debye model of an isotropic solid. It is shown that condensed-matter effects in resonant 
formation are strong, which explains some unexpected experimental results. In particular, the res- 
onance proflles are affected by large zero-point vibrations of the hydrogen-isotope molecules bound 
in the considered crystals, even for high (~ 1 eV) collision energies. This is important for expla- 
nation of the time-of-flight measurements of the //-formation rate, carried out at TRIUMF. The 
calculated mean values of the di/i-formation rate in solid D/T targets, for fixed target temperatures 
and steady-state conditions, are in good agreement with the PSI and RIKEN-RAL experiments. 
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I. INTRODUCTION 



A theoretical study of resonant formation of the muonic molecular ion dtjj, in condensed 
hydrogen-isotope targets is the main subject of this paper. Formation of dt/j, is a key pro- 
cess of muon-catalyzed fusion {/iCF) in a D/T mixture, which attracted particular interest 
because one muon can catalyze more than 100 fusions [1-4] according to the reaction 

dt/i — ^ '^Re + n + ii' + 17.6 MeV . 

Investigation of the /iCF cycle in various hydrogen-isotope targets is also important for 
studies of various phenomena in atomic, molecular, and nuclear physics (see reviews [5-7]). 

Resonant dt/j, formation is due to the presence of the loosely bound state of dt/j, [1] with the 
rotational quantum number J = 1, the vibrational quantum number v = 1, and the binding 
energy ejt,=ii ~ — 0.63 eV. Theoretical methods for calculation of the resonant-formation 
rates were developed for many years (see e.g., Refs. [8-16]). These methods, taking into 
account resonant formation in tfj, collision with one or few molecules, give good agreement 
with the experimental data for dilute gaseous targets. However, such theory is unable to 
explain various phenomena found in experiments with dense fluid and solid hydrogen-isotope 
targets. This concerns a nonlinear dependence of the formation rate on the target density [4, 
17], puzzling temperature effects [18], and the resonance profiles determined by the time- 
of-flight experiments [19-22]. Therefore, it is necessary to consider the influence of many- 
body effects on muonic-molecule formation. In particular, various collective phenomena 
can significantly change this process, which one can expect knowing their role in resonant 
neutron absorption by nuclei bound in condensed matter [23, 24]. 

Condensed-matter effects in resonant neutron absorption can be expressed in terms of 
the single-particle correlation function [24], which has been introduced by Van Hove [25] for 
description of incoherent neutron scattering. This function depends on energy and momen- 
tum transfer to a target and its properties. It is possible to adapt this formalism to the case 
of resonant muonic-molecule formation. 

First estimation of the cZt //-formation rate in solid molecular hydrogens was given by 
Fukushima [26]. He employed a correlation-function formalism, performed ab initio calcu- 
lation of lattice dynamics to determine target properties, and demonstrated an important 
role of phonon processes in resonant dt/j, formation. His calculation was limited to high 
target pressures (~ 10 kbar), where solid hydrogens are classical crystals. However, in /xCF 
experiments, only zero or low pressures (^ 10 kbar) have been applied. As a result, the 
solid-hydrogen targets are quantum crystals with large amplitudes of zero-point vibrations 
of the molecules in the lattices and very different properties. Thus, a special approach is 
necessary to solve lattice dynamics [27, 28]. Owing to this fact and to a rough estimation 
of the transition-matrix elements for dt/j, formation, the results of Ref. [26] are about five 
times greater than the rates determined in the experiments [2, 4]. Moreover, the temper- 
ature dependence of the calculated formation rate, for D2 molecule bound in solid D/T, is 
opposite to what has been recently seen in the RIKEN-RAL experiment [18]. 

A theoretical method of calculating the resonant dd//-formation rate, valid also for low- 
pressure solid hydrogens, has been presented in detail in Rcf . [29] . The correlation function 
used for description of properties of solid polycrystalline D2 has been derived for the Debye 
model of an isotropic solid. The model parameters, such as the Debye temperature and 
the lattice constants, has been taken from the available data including quantum-crystal 
effects [27, 28]. Since the resonances in dd/j, formation on an free D2 molecule are very 
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narrow, their profiles are well-described by the delta function. As a result, the corresponding 
formation rates in a solid are expressed in terms of the same incoherent correlation function 
that is employed for description of incoherent neutron scattering. The theoretical dd/i- 
formation rates lead to the time spectra of dd-fusion products that are in good agreement 
with the data taken at TRIUMF [30]. 

Below we present a method of calculation of the cit/i-resonant-formation rates in con- 
densed hydrogens, for wide intervals of pressure and t/j, collision energy. The profiles of dt/i 
resonances for a free- molecule are described by the Breit-Wigner function [9,31]. In Ref. [24] , 
such profiles have been taken into account for neutron or 7-ray resonant absorption by 
heavy nuclei. It has been assumed that the nuclear mass is not practically changed after 
absorption. As a result, a standard incoherent correlation function was sufficient for de- 
scription of this process. In the case of muonic molecule formation in hydrogens, the mass 
of a target molecule increases greatly after muonic-atom absorption and creation of a small 
muonic-molecular ion. Therefore, we introduce a special resonance correlation function that 
includes this cfi'cct into the target dynamics. Only at lowest collision energies (< 10 meV), 
considered in Refs. [26, 29], an approximation that neglects this mass change can be applied 
since then resonant formation takes place practically in a rigid lattice. Such approach is 
valid for interpretation of experiments performed at lowest temperatures and well-described 
by steady-state kinetics. On the other hand, correct explanation of the time-of-flight exper- 
iments using energetic (~ 1 eV) beams of muonic atoms [20-22] require the knowledge of 
the formation rates at intermediate and higher energies. 

In Sec. II, a brief description of resonant dt/j, formation in an isolated hydrogen-isotope 
molecule is given. A method of calculation of the formation rates in condensed targets, 
using the energy-dependent transition-matrix elements obtained for a single molecule, is 
discussed in Sec. III. In particular, the formulas for the resonant-formation rates in harmonic 
poly crystalline hydrogens are derived. They can be applied to both dt/i and dd/i resonant 
formation. The results of numerical calculations for dt/i formation in low-pressure solid 
hydrogens are presented in Sec. IV. They have been obtained using a full set of the energy- 
dependent transition-matrix elements calculated for the free molecules HD, D2, and DT. The 
dtii formation rates for some typical solid targets are shown as functions of the kinetic 
energy and target temperature. In particular, contributions from different resonances to the 
total formation rates and influence of the ortho-D2 and para-D2 concentration in a target on 
the formation rates are considered. A comparison of the calculated mean rates with some 
experimental results is performed. 



II. RESONANT FORMATION IN A FREE MOLECULE 

First we consider resonant formation of dt/i (the reasoning is analogical for the dd/i case) 
in the following reaction: 

{tfi)F + {BC)i^i,^ [{dtfi)%cee]^^^^ , 
C = H, D, or T and c = p, d, or t , 

where DC is a free molecule in the initial rotational-vibrational state {fiKi) with total 
nuclear spin I. This spin is taken into account for DC=D2. The atom has total spin 
F and center-of-mass (CMS) kinetic energy e. The molecular complex [{dtii)cee] is created 
in the rotational-vibrational state {ufKf) and the molecular ion dt^, which plays the role 
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of a heavy nucleus of the complex, has total spin S. This process takes place due to the 
presence of a loosely bound state of dt/i with rotational number J = 1 and vibrational 
number v = 1. The binding energy released in the reaction above is transferred to 

rotational- vibrational degrees of freedom of the created molecular complex [{dt^)cee]. The 
resonance condition is fulfilled when e takes a specific value This is so-called Vesman's 
mechanism of muonic-molecule formation, introduced in Ref. [8] for the dd^ case. In Fig. 1 
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FIG. 1: Energy diagram for Vesman's mechanism of resonant dt^ formation in collision of a t^, 
atom with an isolated D2 molecule. 



is shown a scheme of energy balance for the )f:/i+D2 case. The formation rate A^^. ^^^^^ 
depends on the elastic width r^^^^^^^x^ of [{dtfi)cee] decay [32, 33] through the channels: 



[{dtfi)%dee]^^j^^ 



> stabilization processes. 



where A/ is the total rate of the stabilization processes, i.e., deexcitations of dtjji and nuclear 
fusion in dtjj,. The value of r^^j^^ ^^^^^^ is given (in atomic units e = H = nie = 1) by the 
equation 

/d^k 
^i^./(^)r5(4-.), (1) 

The transition-matrix element is denoted by Vif{e) and the resonance energy e^jr is defined 
in Ref. [13]. Factor Aif comes from averaging over initial projections and summing over final 
projections of the spins and angular momenta of the system. Vector k is the momentum of 
relative motion of the t/x atom and the molecule DC, connected with kinetic energy e by 
the relation 

e = e/{2M) , 

in which A4 denotes the reduced mass of the system. The general form of Eq. (1) follows 
from the Fano theory of resonant scattering [34] . Integration of this equation over k leads 
to 

A/f 

rSk,..K. = —^A„\V,,(e%)? , kt, = Me",). (2) 
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In the Vesman model, the resonance width is very small, so that the resonant-formation rate 
has the Dirac delta function profile 

^'/K„.fK, = 2nN^,,B,f\V,f{e)\'6{e - e%) . (3) 

where A^moi is the density of hydrogen-isotope molecules in the target. The coefficients Aif 
and Bif in Eqs. (1) and (3) are defined below 



where Wsf — 1 for dt^ and 



W^5F = 3(2F+1)|| I ^1 



(4) 



in the dd^ case. The curly brackets stand here for the Wigner 3j symbol. For asymmetric 
molecules DC, function ^ [Ki) — 1 and in the case of D2 we have 



m) 



I for Ki even 
I for Ki odd . 



A value of factor qd is connected with the number of deuterons in a considered system. 
When dtfi is created in tii collision with an asymmetric molecule DC, = 1, and if D2 
is a target molecule, q^, — 2. In the case of ddii formation in an asymmetric DC system, 
factor qd = 2. For ddfi formation in a D2 target, one has q^ — 4. Coefficient Bif defined 
above differs from that introduced in Rcf. [13] since we omit here the Boltzmann factor 
describing the population of the molecular rotational states in a gas target. We calculate 
the formation rate for a fixed initial rotational state. This rate is however averaged over total 
spin I of the target molecule. If the muonic atoms in the target have a steady kinetic-energy 
distribution f{e,T) at a fixed target temperature T, Eq. (3) can be additionally averaged 
over atomic translational motion, which gives a mean resonant rate A^^^ UfKfC^)- 

Note that, for a given set of the initial and final quantum numbers, the resonance condi- 
tion (cf. Fig. 1) 

e = 4 (5) 
can be fulfilled only when the variable = £11 -|- AE^, is positive. 

In the case of resonant ddfi formation, the rates calculated using Eq. (3) agree very 
well with experiments [35, 36]. On the other hand, the assumption of the delta-function 
profile for dt/i resonances has led to inconsistency with experiments in gaseous D/T targets 
performed at low temperatures [3, 4, 37, 38]. The measured rates are much greater than the 
theoretical predictions based on the Vesman model. It has been pointed by Petrov [9] that 
the dtfi resonances should have broader Breit- Wigner profiles, owing to a finite lifetime of 
the complex. At low temperatures, this leads to significant contributions to the formation 
rates [39] from the subthreshold resonances < 0. Thus, in a general free-molecule case, 
the resonance profile in Eq. (3) can be described by the Breit- Wigner function [9, 11] 

>'u!Ki,UfKf = ^mol Bif I %(£) f _ ^o^p ^ lp2 ' (6) 
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where the total natural width Fs of the resonance is equal to a sum of the effective fusion 
rate A / and the total rate A^^,,^ of back decay of the complex 

rs = Xf + XL ■ (7) 

Equation (6) was employed in Refs. [16, 40] for calculation of dtfi formation rate in a dilute 
D2 gas, which led to agreement with the experimental data [38]. In the limit Fs 0, the 
rate (6) tends to the Vesman form (3). 



III. RESONANT FORMATION IN A CONDENSED TARGET 



A. Method of calculation 



When formation of a muonic molecule takes place in a dense target, it is necessary to 
take into account interactions of the impinging muonic atom with more than one molecule. 
In particular, energy transfer to many molecules is possible, which results in a quasiresonant 
character of the formation process. A quasiresonant mechanism of dtn formation was first 
considered in Ref. [10], for triple collisions t/i+D2+D2, in order to explain a nonlinear 
density dependence of the dtfi formation rate. In this case, formation is possible even if 
the resonance condition (5) is not strictly fulfilled, because an energy excess in the t/i+D2 
system is transferred to the second D2 molecule. The three-body reactions and broadening 
of the resonance profiles were then discussed in Refs. [12, 15, 31, 39]. If a target is condensed, 
it is indispensable to take into account collective motions of target molecules in the process 
of resonant formation. In Fig. 2 is presented a scheme of quasiresonant dtfi formation 




FIG. 2: Energy diagram for quasiresonant formation of dtfi in a D2 molecule bound in a condensed 
target. 



in tjj, collision with a bulk condensed D2 target. The energy balance, including energy 
transfer a; to the target, is shown for the subthreshold resonance corresponding to the 
transition Ui = ^ uj = 2. Since the target molecule and the complex [{dtfi)dee] are bound, 
the corresponding resonance energy Eif is different from the "free" resonant energy 
characterized by the same set of the quantum numbers. 

Owing to a certain analogy between resonant absorption of neutrons and resonant for- 
mation of muonic hydrogen molecules, the methods developed in neutron physics can be 
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adapted for calculation of the rates of resonant ddjj, and dtjj, formation. Resonant neutron 
absorption and emission in condensed targets was first considered by Lamb [23] . His method 
was then generalized by Singwi and Sjolander [24], using the single-particle response func- 
tion Si [25], and applied for description of resonant absorption and emission of 7 ray and 
neutrons in condensed matter. In this Section are derived some expressions for the rate of 
muonic molecule formation in muonic atom collision with molecule DC, bound in a bulk 
hydrogen-isotope target. 

A Hamiltonian J^ot of the system, consisting of a tfi atom in the IS state and a bulk 
condensed DC target, can be written down as follows 

^ot = ^^V| + J^tuiri) + ^DciQi) 

where M^^ is the muonic atom mass and denotes the position of tfi center of mass in the 
coordinate frame connected with the target (see Fig. 3). Operator is the Hamiltonian 



c 




FIG. 3: System of coordinates used for calculation of the formation rate of the complex [(dtfi)cee] 
in a condensed target. 



of a free tfi atom, ri is the tfi internal vector, and .^hc denotes the internal Hamiltonian of 
a free D2 molecule. It is assumed that dt/j, formation takes place in tfi collision with the Ith 
molecule DC. The position of its mass center in the target frame is denoted by R/; is 
a vector connecting the nuclei inside this molecule. Function V stands for the potential of the 
tfi-DC interaction [13] that leads to resonant dtfi formation. Vector Q2 connects the tfi and 
the DC centers of masses. We neglect contributions to potential V from the molecules other 
than the /th molecule because we assume here that distances between different molecules 
in the target are much greater than the DC-molecule size. This assumption is valid for 
condensed hydrogens under low pressure [27, 28]. The kinetic energy e of the impinging 
muonic atom and its momentum k in the target frame are connected by the relation 

e = ky{2M,,) . 

The initial Hamiltonian of the condensed hydrogen-isotope target, corresponding to 
the initial target energy E^, has the following form: 
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where is the position of the jth-molecule CMS (see Fig. 4), Ujj' denotes interaction 
between the jth and j'th molecule, and Mj is the mass of the jth molecule. 




FIG. 4: Position of the impinging tfi atom with respect to the condensed target. 



The coordinate part ip'tot of the initial wave function of the system can be written as 
a product 

iP'tot = ^t^(ri) rDc'iQi) exp(«k ■ Rt,) |0) , (10) 

where |0) stands for the initial wave function of the condensed target, corresponding to the 
total energy Eq. The eigenfunctions of the operators iff^ and i^oc are denoted by ip}^ 

and ipDC^ respectively. Using the relation Rt^ = R; + the wave function ip'tot takes the 
form 

^tot = <^ (ri) rDc'iSi) exp(^k ■ exp(^k ■ R^) |0) , (11) 

which is similar to that used in the case of dt^ formation in a single molecule DC, except 
the factor exp(2k • R^) |0). This factor depends only on the CMS positions of the target 
molecules. 

After formation of the [{dtfi)cee] complex, the total Hamiltonian of the system is well 
approximated by the operator J^^'^^ 

^ot ^ J^tot = '^dti.ir, R) + J^cpUq) + ViQ, r, R) + (12) 

where J^dtfi is the internal Hamiltonian of the dtfi and vectors r and R denote its Jacobi 
coordinates. Relative motion of the dtfi and nucleus c in the complex is described by 
a Hamiltonian ^pix, which depends on the corresponding internal vector g. The final 

Hamiltonian ^ of the target is 



- (13) 



where 

AJ^ = -a -^^Vi , a = 1 - 4^ < - , (14) 



2Mdc ' " Meplx 
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Mcpix is the mass of the complex, and Mqc = Mi is the mass of the DC molecule. A small 
perturbation of potential V, due to replacement of the DC center of mass by that of the 
complex, is neglected here. The eigenfunction and eigenvalue of are denoted by \n) 
and En, respectively. The coordinate part of the final total wave function of the system 
is 

Kt-i^Zir,R)€i^'iQ)\n), (15) 

where ipj^^ and ip^^^xJ stand for eigenfunctions of the corresponding Hamiltonians if^t^ 
and ifcpix- 

For the initial |0) and final \n) target states and for a fixed t/i spin F, the energy- 
dependent (it/x- formation rate ^l$:i,v}Kf{^) calculated using the equation 

^ rs ^ (16) 

{e + E,-e%-Enf + \rf 
with the resonance condition _ 

£ + Eo=4 + F„. (17) 

Equation (16) is analogical to the Breit-Wigner form (6) used for a free molecule. However, 
the transition-matrix element is now given by 

^,0,/n = ('^otl^l'^tot). (18) 

By virtue of Eqs. (11) and (15), the matrix element (18) can be written as a product 

Aiojn = {n\ exp(ik • R0|0) Vif{s) , (19) 

where Vif{e) is the energy-dependent transition-matrix element calculated for a single 
molecule [13]. Averaging the rate (16) over a distribution p„j, of the initial target states 
at a given temperature T and summing over the final target states leads to 

|(n|exp(ik.R0|O)P 



X 



Factor Bif is due to averaging over the initial projections and summation over the final spin 
projections and over the rotational-vibrational quantum numbers. The equation above can 
be written down in the integral form 

^l^KuVfRf = ^mol Bif \Vif\^ Es 

X P^'o |(ri|exp(ik-R,)|0)|^ 

n,no 

dE ^(E-En + Eo) 



^ yS-s^f-Ey + iEl 
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Now we introduce a time variable t to eliminate the 5 function in the equation above and 
then we involve time-dependent operators, which is familiar in scattering theory [41, 42]. 
Using the Fourier expansion of the 6 function one obtains 

1 /"°° ~ 

>^tfK,,u,Kf = ^N^oiB,f\V,j\^rs J dt J2pno\{n\eMi^-'Ri)W^Mit{En-Eo)] 

exp(—iEt) 



n,no 

i-OO 

X / dE 



L 



which, after integration over E, gives 

/oo 
dt exp [-it {e - e^f) - ^Eslt] ] ^ p„„ 

n,no 

X (0| exp(— ik • R;)|n)(n| exp(it£^„) exp(ik • R;) exp(— ii£^o)|0) ■ 

(20) 

The matrix element in Eq. (20) can be expressed as follows 

{n\ exp{itEn) exp(ik ■ R;) exp{—itEo)\0) = (n| exp(it^) exp(ik ■ R/) exp(— 'itJ^)|0) 
= {n\ exp{itJ^) exp{—itJ^) exp(itJ^) exp(ik • R;) exp(— it^)|0) 
= {n\ exp{itJF) exp{-itJ^) exp[ik • R/(i)]|0) , 

where R/(i) denotes the Heisenberg operator 

Ri{t) = exp{itJ^) Ki exp{-itJ^) , 

defined for all / and t. 

Employing the identity |^)(^| = 1 in Eq. (20) we obtain 

^'S<„.,K, = ^nN^oiB.flV.fl' 5,es(k, e - e%) , (21) 

where S^^ is the resonance response function 

5res(k,£-4) = — j dtex^[-it{e-e%)-\rs\t\] 

X ^^res(k,t). 

^es(k, t) denotes here the resonance correlation function defined below 



(22) 



^es(k, t) = { exp[-ik ■ Rz(0)] exp{itJ^) 

exp(-itJf ) exp[ik • Ri(t)])^ , 



X 



(23) 



where (• ■ ■ )t stands for the quantum-mechanical and statistical averaging at temperature T. 

On substitution Jif = and = in the equations above, we recover the well-known 
incoherent response function iSres = Si, which describes incoherent neutron scattering in 
condensed matter [25, 43]. The approximation Jif = Jif is valid when the mass of an 
absorbed particle is much smaller than the mass of a target atom or a molecule. This is 
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a common and good approximation when neutron absorption by a much heavier nucleus is 
considered. However, the difference A^^ between the Hamihonians and cannot be 
neglected in the case of muonic molecule formation since the mass of muonic hydrogen atom 
is comparable with that of hydrogen isotope molecule. 

The partial width -Tf C, ^. k. of back decay of the complex, bound in a condensed target, 
is given by the expression analogical to Eq. (1) 

^!fK,^..K. = SttA/ 1 1^ I Aojnr + E^-s-Eo). (24) 

Using the Fourier expansion of the S function and proceeding as in the case of quasiresonant 
formation process we obtain 

r.fK,,..K. = 27rA/ / ^ |%(.)p5res(k,4 - e) , (25) 

where <Sres denotes function (22) calculated for the initial state |n), with Fs set to zero. 

In order to compare the calculated formation rates with experiments, the summed for- 
mation rates A^. {e) are introduced 

^K,- E ^'5<..,K,, ^. = 0. (26) 

Uf,Kf,S 

In simulations of muon-catalyzed fusion involving energy-dependent rates of various pro- 
cesses the "absolute" formation rates (26) should be used. However, it is convenient to 
consider an effective formation rate A^.(£) that leads to dt fusion in the muonic-molecular 
complex. The fusion probability depends on decay of the created complex into the two initial 
objects: the muonic atom and the hydrogen-isotope molecule. This process competes with 
transitions leading to dt fusion inside the complex. If the lifetime of the complex (< 1 ns) 
is much shorter than its rotational relaxation time, decay takes place back to the initial 
channel. When these times are comparable, it is necessary to include back decay from lower 
rotational states of the complex. In particular, in the limit of very fast rotational relaxation, 
back decay from the ground rotational state Kf = is dominant. Such situation takes 
place in dense targets, where interactions of the complex with neighboring molecules lead 
to fast rotational deexcitation. Calculations presented in Refs. [32, 33] show that rotational 
relaxation of the complex, via scattering on neighboring hydrogenic molecules, is fast at the 
liquid hydrogen density. The effective formation rate is then 

^£ = E ^'i.^fKf , = , (27) 



Kf,S 



where Vf^^ is the fusion fraction 



-PL = ^f/Ts 
and the back-decay rate A^ck is given by 



Abck-E^^^'' ^5^'-E E ^-fKfy.K',- 

F' v[ K[,Kf='a 
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It is assumed here that the vibrational level i^f of the complex is not changed during its 
lifetime. Though calculations of vibrational relaxation of the muonic molecular systems in 
condensed targets have not been performed yet, the available data [28] concerning z/ = 1 — >■ 
relaxation time for H2 in solid (8 /xs) and liquid (12 /xs at 14.2 K) hydrogen suggest that 
such times are much greater than the lifetime of the muonic molecular complex. 



B. Formation in a solid in the strong-binding limit 

Evaluation of the response function 5rcs is difficult, in a general case. The first problem 
is that the operators Rz(t), and AJ^ in Eq. (23) do not commute. However, when 
muonic molecule formation takes place at energies significantly smaller than the mean kinetic 
energy St of molecule DC, the perturbation operator (14) is well approximated by its mean 
value 

AM^ ^ {0\AJ^f^\0) = -a {V%/{2M^c))t ^28) 
= -a (Tt = Asif < . 

Using this approximation in Eq. (23) we obtain 

^es(k, t) ~ CXp {ttAeif){exp[-tk ■ Ri{0)] exp[zk • Ki{t)])^ 
^exp{itAeif)mi{Kt). 

Thus, function ^^cs reduces to the standard incoherent correlation function ^;(k, t) [43], 
multiplied by the factor exp{itAeif ) describing a variation of the mean target energy due to 
its mass change. Hence, the formation rate (21) can be written down as follows: 



/oo 
dt%{k,t) 
-00 



(29) 

X exp [-it{e - Sif) - ^rs\t\] , 



Eif being the resonance energy in the condensed target 



As,f . (30) 



This energy is shifted by Aeif < 0, compared to the free-molecule resonance energy e^j. 
Note that such a resonant-energy shift was neglected in Refs. [23, 24], where absorption of 
neutrons and 7-rays by heavy nuclei were considered. An estimation of the shift in the case 
of 7 emission from a nucleus bound in a solid, similar to Eq. (28) was given in Ref. [44]. 

Using the following relation between ^/(k, i) and the standard single-particle func- 
tion Gs{r,t) [43]: 

^^(r, t) = J d^K exp(-iAt • r) ^ %{k, t) , 

the rate (29) can be expressed as a time and space Fourier transform 

>^'fK„.,K, = N^oiBif \Vif\'J d\dt Gs{r,t) ^^^^ 

X exp [i{K • r — cot) — ^rs\t\\ , 
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where the momentum transfer n and the energy transfer u to the target are 



K = k, 



■if 



(32) 



Analogously, the back-decay width (25) in the strong-binding limit can be expressed 
by Gs{r,t) or by the incoherent response function Si introduced by Van Hove 

Si{K,u!) = — j d^r dt Gs{r, t) exp[i(/« • r — cut)] . 

As a result, Eq. (25) takes a simpler form 

in which 
and 



d'^k 



|y,/(£)|^5,(k,a;'), 



UJ 



'if ^ 1 



(27r)3 

^if = ^if + ^^if ■ 



Ae,f = {n\Aje\n) = - {M^^jM^c - 1) < , 



(33) 

(34) 
(35) 



St being the mean kinetic energy of the complex in the condensed target. 

The equations derived above show that calculation of the formation and back-decay rates 
in the low-energy limit reduces to evaluation of the standard incoherent correlation functions, 
which are well-known in the neutron scattering theory. In particular, for a perfect gas or 
a harmonic solid composed of particles with mass Mmoi, these functions take the simple 
Gaussian shapes [25, 43] 



Gs{r,t) 



M, 



mol 



3/2 



exp 



.2 



_27r7(i) 

For a solid with a cubic Bravais structure, function ^{t) is 

Z{w) 



mol 2 

r 



2M 



mol 



(36) 
(37) 



7(t) 



dw — — ^{coth(i/3rw) [l — cos(wt)] 

— i sin(wt)} , 



w 



(38) 



where the normalized density of vibrational states Z{w) has the following properties: 

/ dw Z{w) — 1 , Z{w) = for W > Wmax , 

Jo 

Z{-w) = Z{w) 



(39) 



and /9t = {k-QT)~^ {k-Q is Boltzmann's constant). 

Solid hydrogens under low pressure, used for studies of muonic atoms and molecules, are 
quantum molecular crystals. They have the Bravais fee polycrystalline structure or the hep 
polycrystaUine structure [27, 28], for which Eqs. (36) and (37) are fair approximations. As 
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a result, on substitution M^oi — Mdc, we obtain the phonon expansion for the resonant- 
formation rate 



SF 
ViK. 



N^oiBif\Vif\^eM-'2W) 



X 



4 o 



(21^)" 



ni 



(40) 



in which 



9ri{w) 



grn[w) 



1 

2^ 



dz ■ 



2 I 1 7-'2 



gi{z + w,T) , 



dw' gri{w - w') gn-i{w') , 



(41) 



and 



Z{w) 



7(00) 



w 



[^bW + i] , 



(42) 



/oo 
dw' gi{w - w') g„^i{w') , 
-co 

/oo 
dwgn{w) = 1 . 
-CX) 

The exponent 2W of the Debye- Waller factor exp{—2W), famihar in the theory of neutron 
scattering, is 

^2 

2W{k'^) = 7771—7(00) 



°° , Z(w) ^ , 

aw — — - coth (^prwj , 



2Mmol 7o ^ 

where 7(00) stands for the limit of ^{t) at i — > 00. Function rij^iw) denotes the Bose factor 

- [(^MPtw) - (43) 

The Breit-Wigner term in expansion (40) describes recoil-less resonant formation. The 
sum with higher powers of 2W correspond to quasircsonant muonic molecule formation with 
simultaneous phonon creation or annihilation. In particular, the term with n = 1 describes 
formation connected with creation or annihilation of one phonon. In the strong-binding limit 
21^ <^ 1, only few lowest terms in expansion (40) are significant. The phonon expansion (40) 
is more general than an analogous expansion in Ref. [24], which includes the Breit-Wigner 
factor only in the nonphonon term. This factor should be taken into account also in the 
phonon terms, unless the natural resonance width is much smaller than Wmax- For 2W > 1, 
the approximation (28) and Eq. (40) are no longer valid. 

When dd/i formation is concerned, the resonances are very narrow. Thus, in this case, 
the limit — > is practically reached. The Breit-Wigner factor tends to the 5-function 
profile and gm 9n- As a result, Eq. (40) takes for dd/i a simpler form, derived in Ref. [29] 



X 



{2wy 



n=l 



nl 



(44) 
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A phonon expansion can also be applied for estimation of the back-decay rate. After 
integration of Eq. (33) over direction of k one obtains 



A, 



dk k^\Vif{e)\^Si{k\uj'). 



Substitution of the phonon expansion for Si into equation above and then integration of the 
5-function term lead to 



SF' 



UfKf,UiKi 



A. 



TT 



OO 

+ V / dkk'^\Vif{e)\'^exp{-2W) 

n=l 



(45) 



in which 



2W = 



2M 



cplx 



■7(00), 2Wif^2W{hf), 



and kif — ^J~2Mi^f 



(46) 



are calculated for the harmonic lattice with the bound muonic molecular complex. Note 
that Eq. (45) is valid only if a main contribution to the integral comes from small k. 



C. Formation in the weak-binding limit 

When the incident momentum of the muonic atom is large, the formation time of a muonic 
molecule is short compared to the characteristic time scale of the dynamic response of the 
bulk target. Thus, a contribution to the response function (22) from short times is dominant. 
As a result, it is sufficient to keep only linear terms in t while evaluating an asymptotic form 
of the correlation function ^es(k, t). In calculations, we shall use the following operator 
relation: 

exp(i) exp(S) = exp(i + S + (7) , (47) 

where 

C = 5] + ^ [[i, 5], 5] + ^ [[^, i], i] 



+ J- 

~ 24 



+ 



Operator (7 = only if A and B are commuting operators. 

The operators AM' and J^, defined by Eqs. (9) and (13), do not commute and the 
operator C in the expression 

exp{ii(jr + AM')} ex.p{-itJf) = exp{itAje + C) 
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turns out to be a sum containing higher powers of t. Since in this approximation we restrict 
to terms hnear with respect to t and to the parameter a; < |, the operator C in the relation 
above can be neglected and thus the correlation function takes the form 

^resiK t) = ( exp{-^k ■ Ri{0)} cxpiitAM') 
xexp{ik-Ri(i)}>^. 

Now we involve the basic approximation 

RKt)^R(0) + (P;/MDC7)t, (49) 

where P; denotes the momentum operator of the Ith molecule. This approximation is valid 
for t — > 0. After substitution of Eq. (49) in Eq. (48) and multiple use of the Eq. (47) we 
have 

^es(k, t) ~ exp (it — - — I / exp ( —ita^^ — 
X ( exp tt— — 

Since the argument of the second exponential is small, we can use the following approxima- 
tion: 

f ■ Pi W f ■ / Pi \ 

exp —ita — ) a; exp —ita ( — ) 

— ex.p{itAeif) 

which involves the resonance-energy shift (28). Substitution of the above equations 
in Eq. (22), with the definitions (30) and (32) taken into account, leads to 



Sres{i^,u}) = ^ / dt exp 



—iujt — ■krs\t\ -\- it 

■cplx. 



X ^exp ^it 



2 

/*Pz 



2M„ 



M, 



cplx 



(50) 



When the motion of the molecule DC is well described by an isotropic harmonic potential, 
the Bloch identity 

(expg)T = exp(|(Q2)^) (51) 

may be applied for an operator Q being a linear combination of the Bose operators of 
creation and annihilation. Since momentum can be expressed by such operators (see 
e.g., Ref [43]), we have 

In the case of cubic symmetry. 



^c'plx 
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and this is a fair approximation even for other lattices. Thus 



8 M, 



DC 



/ Pf \ _jI 



3 Mepbc \2MDc/r2Mepix' 

which finally gives the following Doppler width: 



I — — St ojr 1 



3 M 



cplx 



(53) 



with the recoil energy 

a;R = KV(2Mepix). (54) 
In the case of a solid hydrogen target, the mean kinetic energy of the bound molecule equals 



St 



dw Z{w) w [rZg (w) + |] . 



(55) 



This energy is much higher than St — ^kBT for a corresponding Maxwellian gas, unless the 
temperature is sufficiently high. This phenomenon was first taken into account by Lamb [23] , 
in resonant neutron absorption in solid crystals. In particular, for a low-pressure solid or 
liquid deuterium, St ^ 5 meV [45] due to a large zero-point motion of D2 molecules in 
a given target. The effective target temperature Tes corresponding to St is then defined as 



For the considered solid D2 case, Tes ~ 40 K. 

Substitution of Eqs. (52) and (54) in Eq. (50) leads to 



(56) 



S^es{i^, ^) ^ TT (it exp [ - i(a; - ujR)t 
27r J_^ 



(57) 



krs\t\ 



4^rcs'' J • 



Then, applying the convolution theorem to the Fourier transform of a product, we obtain 
the asymptotic form of the resonance response function 



5res(K,a;) 



Ps 



dz 



X exp 



1 r2 

4^ s 



Z + UJ — 



(58) 



By virtue of Eq. (58), the formation rate (21) in the weak- binding limit takes the form 



A 



SF _ AT n IT/. |2 -^'S' 

.,K.,.fKf-^'^mol-tlif\Vif\ ^^^^ 



X 



dz 



^2 I 1 7-'2 



exp 



Z + UJ — UJ-R^ 



(59) 
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This equation is similar (apart from the muonic- molecule factor iVmoi-Bj/l^/P) to the expres- 
sion for resonant absorption of neutrons in a gas target, obtained by Bethe and Placzek [46]. 
However, the resonance width (53) and recoil energy (54) take into account a change of the 
target particle mass in the absorption process, which is neglected in their work. 
In the limit Fs — > 0, Eqs. (58) and (59) tend to the following expressions: 



TT 



exp 



^res J 



(60) 



and 



ifl 



X 



exp 



(61) 



respectively. Function (60) has the Gaussian form, identical with that used for description 
of incoherent scattering at large energies. However, the Doppler width (53) and the recoil 
energy (54) in 5res are different from the corresponding variables 



(62) 



and 



CUR = kV(2M^oi) , (63) 
which determine the asymptotic form of the standard incoherent response function Si [43]. 
Function S^es tends to Si if the approximation Mcpix ~ Mdc is valid. However, in the case 
of muonic molecule formation, this is only a rough approach because the mass of a muonic 
hydrogen atom is comparable with the mass of a hydrogen isotope molecule. Note that, in 
the strong-binding limit, mass Mcpix enters only the resonance-energy shift (28). The phonon 

expansion (40) is expressed in terms of ~ k^/2Mj)c, not in terms of ~ K^/2Mcp\x. 
The reason is that Eq. (40) is valid for small collision energies, when the target molecule is 
strongly bound in the lattice. Therefore, the momentum is mostly transfered to the whole 
crystal. 

Function (60) can be used for evaluation of the back-decay rate, if large final momenta 
give main contribution to the integral (25). After integration over direction of k in Eq. (25), 
with the asymptotic function (60) inserted, one obtains 



r 



SF' 



Ai 



if 



VfKf,ViKi 



dkk'\Vif{e 



X exp 



UJ 



LU 



R 



(64) 



where uj' is defined by Eq. (34). The parameters A-^^s f^^nd uj'i:, are calculated from Eqs. (53) 
and (54), using the replacements A'/dc ^ Afcpix and St <^t- 

Let us note that Eqs. (58) and (60) are general since they are derived in the impulse 
approximation (49) without using specific properties of a given target, apart from the single 
parameter St- Therefore, they are valid for liquid and dense gaseous hydrogens. They can 
be also used for description of resonant absorption processes other than muonic molecule 
formation, when a mass change cannot be neglected. In such a case, Fs should be replaced 
by an appropriate natural resonance width. 
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D. Formation in a solid at intermediate energies 



The formation rate calculated according the asymptotic form (59) becomes very inaccu- 
rate when the collision energy is much greater than the maximal frequency Wmax of a crys- 
tal. In particular, this concerns recoil-less formation, which is dominant at lowest energies. 
Therefore, at the intermediate energies, it is reasonable to represent the formation rate as 
a sum of the exact nonphonon term from expansion (40) and the subsequent phonon terms, 
which we obtain below in the impulse approximation. Using Eqs. (47) and (51), it can be 
shown that the relation 



X exp 



\a[it- l{a + 2)^Tt'^] 



2M, 



DC 



(65) 



is valid in this approximation. Inserting Eqs. (37) and (65) into Eq. (22) we have 

1 



2n 



dt exp< —ituj — iFslt 
I 

+ [--f{t)+iat 

-la{a + 2)^Tt'] 2^ 



(66) 



K 



DC 



Substituting the approximation j(t) ^ — it + ^ Srt'^ for short times into Eq. (66) and 
integrating over t yields the asymptotic form (58) of the response function. However, we 
now expand Eq. (66) in powers of 



1 ^ (2wy 

S,,,{K,uj)= — exp{-2W)J2 

n=0 

f oo 



X 



/oo 
dt exp{-iut - ^rs\t\) [^{t)Y, 
-oo 



(67) 



^{t) = 1 + i 



a 



a 



7(00) 



2 (1 

3 7(cxo 



t^ ■ 



Function Qn is defined by Eq. (42). The integral over t is estimated using the exponential 
approximation to function ^ 



^{t) exp(x) , 
where x contains only leading terms in t and 

■7(00) , 



it 

la 



T- ~ ""' _ 1 /\2 /2 



(68) 



Mepix 

la = 



Ml 



DC 



^2^4 Mdc ^ 

" " 3 Mcplx la 



7^ 

la 



Then, integration in Eq. (67), using the convolution theorem, leads to 



5res(«;,Cu) = exp(-2iy) 



27r + |r| 



{2wy 



(69) 



n=l 
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where 



1 f°° 



X [^^(z + a;^) + 1] , 



1 + a ' 



and, for n > 2, 



^2 I 1 n2 



X exp 



{z + u ~ nj'^af 
2nAl 



The first term of Eq. (67) has been replaced in (69) by the exact Breit-Wigner term. Also the 

one-phonon (n = 1) contribution to iSrcs, is replaced here by a more accurate term depending 
on 01. Function gi is calculated on substitution of the exact function 7(t) for a harmonic 
solid into Eq. (66). Every multiphonon term in Eq. (69) is represented by the convolution 
of the Breit-Wigner profile with a Gaussian obtained using Eq. (68). It now follows that 



rs 

n=l 



X 



+ 27r^0„(a;) 



{2wy 



(70) 



The form of this expansion is similar to that of Eq. (40), derived in the strong-binding hmit. 
However, functions 0„ are obtained in the impulse approximation and they are different 
from the corresponding functions gm given by Eq. (41). For the one-phonon term, we have 
01 (cj) = grii^^a), which is the direct result of using the exact 7(t) in derivation of 0i. Thus, 
Eqs. (70) and (40) give the same rate at smallest energy transfers. At large s, when many 
multiphonon terms are important, the target response no longer displays a rich structure. 
The rate (70) tends therefore to the simpler form (59), which is characterized by the recoil 
energy (54) with the correct mass Mcpix. 

In the limit Fs — >■ 0, the rate (70) takes the form similar to Eq. (44) 



X 



{2wy 



n=l 



(71) 



with the expansion coefficients 



0i(u;) = ^K(u;«) + l], 



0n(<^) = 



(27rn)V2z\^ 



exp 



[UJ 



2nAl 



n>2. 



The back-decay rate can be calculated analogously. The result is given by Eq. (45) with 
functions Qn replaced by the corresponding functions Qn from Eq. (71). It is also necessary 
to make the following substitutions: Mj^c ^ -^cpix and St — > St- 
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IV. RESULTS OF CALCULATIONS FOR dt/j, FORMATION IN SOLID HYDRO- 
GENS 



In this Section, the rates of resonant dt/j, formation in sohd HD, D2, and DT are calculated. 
It is assumed that these targets are kept at zero or low pressures (<S 10 kbar), which 
corresponds to the TRIUMF or RIKEN-RAL experimental conditions. Measurements of 
the formation rates at TRIUMF have been performed using energetic (~ 1 cV) beams of 
t/j, atoms. Therefore, the rates are evaluated here in a wide energy interval £ ^ 1 eV. This 
involves resonant dt/i formation with simultaneous excitations of a few lowest vibrational 
levels of the muonic-molccular complex. The values of the rates are given for a normalized 
target density of 4.25 x 10^^ atoms/cm"^ (liquid- hydrogen density). 

Since exact forms of the vibrational-state distribution Z{w) for the experimental poly- 
crystalline targets are not known, the Debye model of an isotropic solid has been used in 
the calculations presented below. The values of the Debye temperature ©d are taken from 
the available literature [27, 28]. 

The resonance energies and energy-dependent transition-matrix elements for isolated tar- 
get molecules HD, D2, and DT, calculated according to the method presented in Ref. [47], are 
the starting point for evaluation of the formation rates in solid hydrogens. The transition- 
matrix elements are available for the rotational transitions Ki = 0,1 ^ Kf = 0, . . . , 9. 

Resonant dt^ formation in a bound D2 molecule is the most complicated case. The lowest 
resonances, corresponding to the vibrational transition i/j = — > i// = 2 and different rota- 
tional states Ki and Kf, are located in the vicinity of £ = with the radius of a few tens meV. 
The resonance energies in this region, for a free D2 molecule and for a D2 bound in a 3-K 
solid deuterium, are shown in Table I. In particular, there are several subthreshold reso- 



TABLE I: Resonance energies for dt^ formation in t^x scattering from a free D2 molecule {s^j) and 
from a 3-K solid-D2 target (ei/), corresponding to the vibrational transition i^j = — > i^/ = 2. 
These energies are given in the corresponding center-of-mass systems. 
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nances that give significant contributions to the low-energy rates, because of wide resonance 
profiles. The resonance-energy shift (28) for a deuterium target at 3 K is ^Sif = —2.29 meV. 
Resonances in the upper spin state F = 1 have much smaller energies than those for F = 
with the same rotational quantum numbers. In particular, the largest values of Sif for F = 1, 
shown in Table I, are due to the excitations Ki = 0, 1 ^ Kf = 4. The only matrix elements 
that do not tend to zero at e ^ correspond to the dipole transitions Ki = Kf = 1 and 
Ki = 1 ^ Kf = 0, 2. For F = 1, all these transitions are associated with Sif < —50 meV 
and thus they give very small contribution to the resonant (it/i- formation rate. As a result, 
the low-energy rate is determined mainly by tfi scattering in the F = state. However, even 
for F = 0, the dipole transitions are connected with negative resonance energies, though 
much closer to e = than in the F = 1 case. The lowest positive resonances appear in 
the transitions Ki = Kf = 3, 4 and Ki = 1 —>■ Kf = 3, 4. They are characterized by 
strongly varying transition-matrix elements [47], which is illustrated in Figs. 5 and 6. Let 




tjo, energy (meV) 

FIG. 5: Matrix elements |l^/(e)p versus tfx energy for the transitions Ki = ^ Kf = 0, 1, 2, 3, 4 
and = — > z^/ = 2. The vertical lines denote energies e,/ of the lowest resonances. Labels 
"i — > /" stand for the rotational transitions Ki Kf. 
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FIG. 6: Matrix elements versus t/i energy for the transitionsATj = 1 — > Kf = 0, 1, 2, 3, 4 

and Vi = Q ^ Vf = 2. Notation is the same as in Fig. 5. 
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us note that this situation is very different from the ddn case, where low-energy formation is 
determined by the dipole transitions, with the matrix elements slowly varying below a few 
tens meV [29, 47]. Another difference between the ddfi and dtfi case is involved by larger 
separations of the neighboring dt^ resonances corresponding to i^j = and Ki = 1. There- 
fore, for dt^ one can expect more pronounced differences between resonant formation in solid 
ortho-D2 and para-D2 than those found for dd^ case [48] . Most pure-deuterium experiments 
in /iCF have been carried out in targets with the statistical mixture of ortho and parastates 
(called "normal" deuterium nD2, according to the nomenclature used in Ref. [28]). 




FIG. 7: Transition-matrix element |VJ,j(e)p for resonant dtfj. formation (transition = — > z^j = 2, 
j = 1 — > Kf = 3, dashed line) and the response function S^esii^, £ — ^if) (in arbitrary units, solid 
line) for the resonance ^ = 0^5 = 1 in 3-K para-D2. The peak of the Breit-Wigner term from 
Eq. (69) is centered at the resonance energy Sif = 2.7 meV. 



In Fig. 7 is shown |V^,/(£)P for the transition Ui = Uf = 2, Ki = 1 Kf = 3, together 
with the response function (69) for the resonance F = —>■ S = 1 located at = 2.7 meV. 
The phonon terms in iSres are calculated assuming = 0, since in this example we want to 
neglect their convolution with the Breit-Wigner profile. There is a strong contrast between 
resonant formation of the molecules dtfi and ddfi [29] in a solid deuterium. In the dt/j, 
case, the wide Breit-Wigner peak is not so much pronounced as the narrow recoil-less dd/j, 
resonances. The matrix element |Vi/(£:)p raises by a few orders of magnitude within the 
width of 100 meV of the multiphonon distribution. Thus, the phonon contribution to the 
(iiyLi-formation rate is comparable with the nonphonon one, already above a few meV. This 
means that a detailed form of the density Z{w) of vibrational lattice states is necessary for 
accurate calculation of the low-energy dt/i-formation rate in a solid D2. A shape of the 
phonon spectrum in the energy-dependent rate is strongly distorted, which one sees in Fig. 8 
evaluated using Eq. (70). Nevertheless, the one-phonon and two-phonon terms are clearly 
distinguished in the curve corresponding to para-D2. In ortho-D2, the resonance with the 
lowest £j/ > is located at 10 meV. Therefore, the Breit-Wigner peak is strongly suppressed 
by the Debye- Waller factor and the rate is quite fiat. At e — > 0, the rates are determined 
by the wings of the Breit-Wigner peaks, because phonon contribution to the rates vanishes 
when K approaches zero. For F = 1, the main resonances are far from the considered 
low-energy interval (see Table I). The rates shown in Fig. 9 are thus determined by the 
Breit-Wigner wings of the deep subthreshold resonances with small contributions from the 
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FIG. 8: Low-energy dt/i- formation rate for F = in a 3-K solid nD2 ("normal" deuterium [28]), 
ortlio-D2, and para-D2, calculated using Eq. (70). 
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FIG. 9: Low-energy (it^-formation rate for F = 1 in the same targets as in Fig. 8. 



weak resonance {Kf = 5) located a.t Sif > 0. As a result, the formation rates for F = 1 are 
lower by two orders of magnitude than those for F = 0. 

Resonances in tfi scattering from D2, corresponding to the vibrational excitations z/j > 3 
of the [{dtfi)dee] complex, are located at higher energies e > 0.2 eV. Therefore, they are well 
described by the asymptotic form (59), which is independent of Z{w). Therefore, the forma- 
tion rate is determined accurately using only the mean kinetic energy of a D2 molecule. 
The formation rate in a 3-K solid nD2 is plotted in Fig. 10, for several Vf. For comparison, 
in Fig. 11 is shown the (it/i- formation rate for 3-K gaseous nD2. The energy-dependent rate 
for a perfect deuterium gas has been calculated assuming a 3-K Maxwellian distribution of 
the D2 kinetic energy. This rate includes only formation due to two-body t/i+D2 collisions. 
The resonant-formation rates presented in Figs. 10 and 11 display a striking difference be- 
tween the gas and the solid case. At e — > 0, the theory developed for two-body collisions 
in a perfect gas gives a negligible resonant formation rate. This result disagrees with the 
average formation rates determined by measurements performed in liquid and cold dense-gas 
targets [3] and in solid [18] targets. The rate for the solid shows a strong contribution from 
the subthreshold resonances, which leads to a large rate in the limit e ^ 0. Solid-state 
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FIG. 10: Resonant dt/i-formation rate in a 3-K solid nD2 for F = and F = 1. The label I'f 
denotes the vibrational state of the created [{dtfi)dee] complex. 
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FIG. 11: Resonant dt/i- formation rate in 3-K gaseous nD2 calculated in the laboratory frame. 



effects are also significant at higher energies. The resonance peaks in the solid are much 
broader than those in the gas because of a large effective target temperature connected with 
strong zero-point motion of D2 in solid deuterium [45]. The widths of the peaks increase 
with rising recoil energy. However, the centers of higher-energy peaks in the both targets 
have similar locations since, in the impulse-approximation limit, the recoil energy (54) for 
the muonic-molecular complex bound in a solid equals to that for the isolated complex. 
A small difference of the resonance energy between the solid and the gas is negligible 
for £4/^1 meV. 

Calculation of the (it/i-formation rate for a solid HD or DT is simpler than for D2 since in 
the HD or DT case there are no significant resonances in the close vicinity of e = 0. This is 
caused by different values of the rotational and vibrational quanta for these three molecules. 
The HD molecule is the lightest one and the resonances connected with uj = 2 are situated 
in HD above 0.1 eV [14, 47]. As a result, contributions from various multiphonon processes 
to the formation rate plotted in Fig. 12 cannot be distinguished. The resonance peak for 
Uf = 2 in HD is the strongest dtfi resonance found for the three considered molecules. 
In Fig. 13, the diyLi-formation rate is shown for a 3-K solid DT target. The lowest peaks. 
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FIG. 12: Resonant /Li- formation rate in 3-K solid HD. 
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FIG. 13: Resonant dt/i-formation rate in 3-K solid HT. 




which aheady take the asymptotic form (59), correspond here to Uf = 3. The rotational and 
vibrational quanta are smallest for DT, so that the main (lowest Kf) resonances connected 
with Uf = 2 are located deeply below e = 0. Thus, a contribution to the formation rate from 
the subthreshold resonances is very small and is not apparent in this figure. At 3 K, the 
effective target temperature Tgfj, determined by Eqs. (56) and (55), equals about 41 K for HD 
and 50 K for DT. The resonance shift ^^j/ obtained from Eq. (28), equals —2.71 meV in 
the case of HD and —1.97 meV for DT. 

The dtfi resonances in solid HD and D2 were directly observed at TRIUMF [19-22] using 
the energetic t/x-atom beam and time-of-fiight techniques. However, Monte Carlo simula- 
tions (see e.g., Ref. [49]) were employed for interpretation of the experimental data. Such 
a procedure was indispensable since the time-of-fiight spectra cannot be uniquely inverted 
because of the geometry used and the energy loss of t/x atoms in the reaction layer, prior to 
resonant formation of the muonic-molecular complex. In those simulations, the calculated 
c?t/i-formation rates for 3-K gas targets (such as that shown in Fig. 11) were applied because 
the theoretical formation rates for a low-pressure solid were not available. A detailed anal- 
ysis of the data was performed by Fujiwara [19]. He found more (it-fusion events at lowest 
and highest t/i energies than it was predicted using the perfect gas model. Much broader 
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resonance peaks, which we present in Fig. 10, can certainly improve the fits to the TRIUMF 
data. Also, the analysis of the fusion-product yield [19] proved that the low-energy dt^- 
formation rate in solid deuterium was much higher than that predicted by the two-collision 
gas model. In particular, this concerns formation for the state F = 1. The theoretical rates 
presented in Fig. 9 support this finding. 

A two-peak structure of the calculated time-of-flight spectra for dtfi resonances in HD, 
obtained assuming a 3-K gas model, was not confirmed by the solid-HD data [21]. How- 
ever, one may expect much better agreement with the TRIUMF data when the rate shown 
in Fig. 12 is used instead of the very pronounced peaks evaluated for a 3-K HD gas. A possi- 
bility of wider resonance peaks with a fixed Doppler width of 50 meV was already considered 
in Ref. [21], which did not give good fits to the data. Such a result is now clear since, accord- 
ing to Eq. (53), the Doppler width of a resonance in a condensed target increases with the 
rising recoil energy cur. Simultaneously, the resonance height (59) decreases for higher 
so that the resonance strength is preserved. 
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FIG. 14: Rate of resonant dtfi formation in tfi(F = 0) scattering from a D2 molecule bound in 5-K 
and 16-K solid D/T(Ct = 0.4) targets. 
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FIG. 15: Rate of resonant dt^ formation in tii{F = 0) scattering from a DT molecule bound in 
5-K and 16-K solid B/T{Ct = 0.4) targets. 
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In Figs. 14 and 15 are shown the resonant-formation rates for the molecules D2 and DT 
bound in a solid D/T target. An equilibrated mixture of the molecules D2, DT, and T2 is 
assumed, for the tritium isotopic concentration Ct = 0.4. Target temperatures 5-16 K were 
applied in the RIKEN-RAL experiment, in which an unexpected temperature dependence 
of the (it/i-formation rate in solid D/T mixtures [18] was found. The corresponding target 
density is almost constant. A similar hydrogens mixture, kept at 15 K, was also used in the 
PSI experiment [4]. In the both experiments, time spectra of neutrons from dt fusion were 
measured. The data were interpreted using a standard steady-state kinetics, assuming that 
tfi atoms were thermalized. Formation from the state F = 1 is negligible for an appreciable 
tritium concentration as the spin-flip transition F = 1 ^ in low-energy tfi + t collision 
is very fast [50]. The theoretical energy-dependent rft/i-formation rates display a weak 
temperature dependence. One can expect such a behavior since, for any temperature of 
a low-pressure hydrogen-isotope solid, the limit T/O-o -C 1 is achieved {O-o ~ 100 K) 
and changes of 0^ are very small [27, 28]. As a result, the response function (69) and 
thus the formation rate (70) are always close to their limits for T/^d — > 0. Therefore, 
changes of the average formation rate, determined using steady-state conditions, can only 
be ascribed to different t/i-energy distributions corresponding to various target temperatures. 
An accurate comparison of the theory with data requires Monte Carlo simulations of the /xCF 
cycle in a given solid D/T mixture, which can be performed in future after completion of 
a full set of the differential cross section for muonic atom scattering in mixed D/T crystals. 
The t/i-energy distribution in steady-state conditions is a crucial information. A shape 
of such a distribution is non-Maxwellian and the mean t/x energy is greater than ^k-oT, 
due to solid-state effects and a possible admixture of epithermal t/i's from the reaction 
dfx + t ^ + d and from back decay of the muonic-molecular complex. The latter effect 
was studied in Refs. [38, 51] with the use of Monte Carlo simulations, in the case of gas and 
liquid targets. In a high-density target with medium or high Ct, this effect is small, which 
is confirmed by the PSI fits [4]. 

Averaging the energy- dependent rate from Fig. 14 over the t/x-energy distribution leads 
to the mean resonant rate shown in Fig. 16. The energy distribution of tfi atoms, being 
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FIG. 16: Mean rate of resonant dtfi formation in tfi{F = 0) scattering from nD2 molecules bound 
in solid D/T (Ct = 0.4) as a function of the target temperature. The dashed line represents the 
same rate scaled by the factor Sx = 0.86. Also is shown the result of PSI measurement [4] for 
a similar target (T = 13 K, (/? = 1.45). 
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in thermal equilibrium with phonons, is assumed to be proportional to Z{e) nB{e,T). The 
average tfi energy obtained using this function ranges from 1.2 meV for T = 5 K to 3.4 meV 
for T = 16 K. It is evident that the rise of the formation rate, above about 3 K, is mainly due 
to t/j, entering into the region of the recoil-less resonant peak in para-D2, centered at 2.7 meV. 
Phonon processes in both ortho-D2 and para-D2 lead to a smaller rise of the rate. The 
calculated formation rate is close to the PSl result for T = 13 K [4]. A coincidence of the 
theoretical curve with the data is obtained, as in the case of the TRIUMF measurements [20] , 
on scahng by the factor 5'a < 1, which can be ascribed to inaccuracy of the calculated 
transition-matrix elements. Here, we have Sx — 0.86, which is consistent with the result 
of Ref. [20]. 

In the RIKEN-RAL experiment [18], about 20% decrease of the fiCF effectiveness has 
been found for the target-temperature change from 16 to 5 K, independently of the tritium 
concentration. In order to explain this effect, several hypotheses have been considered. 
The hypothesis of a significant change of the mean resonant dtn formation rate A°^^ (for 
F = 0) has led to best fits to the data. Kawamura et al assume that the two components 
of A[]t^, namely the rate A[)^|^^ of resonant formation for the D2 molecule and the analogous 
rate A|]^^^ for the DT molecule, are comparable. At 16 K, they use A|]^^^ = 3.5 x 10® s~^ and 
~X^fJ = 1.6 X 10® s-i [18]. All temperature dependence of A^,^ = Cd A°;°'+a ~X^fJ {Cd is the 
deuterium isotopic concentration) is ascribed only to X^dt^^- Other rates in the steady-state 
kinetics being fixed, about 30% decrease of Xdt^^ between 16 K and 5 K has been obtained. 
Thus, for Cf = 0.4, the respective change of A^^^ equals about 25%. This finding agrees quite 
well with analogous 20% decrease of the theoretical rate plotted in Fig. 16. However, theory 
predicts that the low-energy rate A°J^^ should be smaller by a few orders of magnitude than 

the corresponding rate A°J^^, since the strong resonances in tyU+DT scattering are far from 
the region e ^ 0. Averaging the rate presented in Fig. 15 over the t/x-energy distribution 
gives A°J^ = 2.6 x 10^ s^^ This value agrees well with the rate A°;J^^ = (1.8±0.7) x 10^ s-^ 
determined for a 30-K liquid D/T in the PSI experiment [4]. Note that the formation rate 
in the solid is somewhat greater than the corresponding rate in the liquid, which is a general 
law confirmed by experiments. Thus, according to the presented calculation and to the PSI 
results, A|]^^ ^ ^^t^^ ■ This means that in the steady-state analysis of Ref. [18], a somewhat 

greater value of A|]^^^ should have been assumed. In fact, Monte-Carlo simulations similar 
to that performed for gaseous D/T [38] are indispensable for an accurate analysis of such 
experiments, since several rates change significantly at lowest energies and thcrmalization 
process of muonic atoms in solid hydrogens is complicated [29, 52, 53]. It depends on the 
target temperature, isotopic concentration, and rotational population. A full set of the 
differential cross sections for muonic atom scattering in mixed sohd D/T is necessary for 
accurate description of /iCF in such a target. 

V. CONCLUSIONS 

A method of calculating the rates of muonic-molecule resonant formation in collision of 
muonic atoms with condensed hydrogens has been developed. In the case of polycrystalline 
hydrogen-isotope targets, detailed calculations have been performed using the Debye model 
of an isotropic harmonic solid. Vahics of the resonant- formation rates have been computed 
for resonant dt/i formation in frozen D/T and HD targets, for collision energies < 1 eV. 
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These rates arc very different from those obtained for dilute gaseous hydrogens and exhibit 
strong soUd-state effects. 

At lowest energies, contributions to the total rate from formation in a rigid lattice and 
from formation with simultaneous phonon processes can be distinguished. In the high- 
energy region (e ^ 0.1 cV). for any target, the rate takes a general asymptotic form, which 
depends on the mean kinetic energy of a target molecule. For low-pressure solid and liquid 
hydrogens, this energy is much greater than the corresponding energy in a perfect gas. As 
a result, condensed-matter effects in resonant formation do not disappear even at highest 
collision energies. Since the main dt/j, resonances for HD and DT are located far from zero 
energy, in these cases it is sufficient to use only the asymptotic expression for the formation 
rate. 

The calculated resonance profiles in solid are much broader than in the dilute-gas case. 
Experimental evidence supporting this conclusion has been found in the time-of-flight mea- 
surements of dtfi resonances at TRIUMF. A quantitative comparison of the theory with 
these experiments requires however complicated Montc-Carlo simulations. 

The mean values of the dt/x- formation rates for D2 bound in the sohd D/T mixtures, 
averaged over the t/j, kinetic energy under the steady-state conditions, agree well with the 
PSI and RIKEN-RAL data. Also a temperature dependence of the mean formation rate, 
determined at RIKEN-RAL for temperatures 5-16 K, is revealed by the theory. 



Acknowledgments 

We would like to thank Profs. K. Nagamine and L. I. Ponomarev for supporting this 
research and Drs. M. C. Pujiwara and G. M. Marshall for valuable discussions. 



[1] S. S. Gershtein and L. I. Ponomarev, Phys. Lett. B 72, 80 (1977). 
[2] S. E. Jones, A. N. Anderson, A. J. Caffrey, et al., Phys. Rev. Lett. 56, 588 (1986). 
[3] W. H. Breunlich, M. Cargnelli, P. Kammel, et al., Phys. Rev. Lett. 58, 329 (1987). 
[4] P. Ackerbauer, J. Werner, W. H. Breunlich, et al., Nucl. Phys. A 652, 311 (1999). 
[5] W. H. BreunHch, P. Kammel, J. S. Cohen, and M. Leon, Ann. Rev. Nucl. Part. Sci. 39, 311 
(1989). 

[6] L. I. Ponomarev, Contemp. Phys. 31, 219 (1990). 
[7] P. Proelich, Adv. Phys. 41, 405 (1992). 

[8] E. A. Vesman, Pis'ma Zh. Eksp. Teor. Fiz. 5, 113 (1967), [JETP Letters 5, 91 (1967)]. 

[9] Yu. V. Petrov, Phys. Lett. B 163, 28 (1985). 
[10] L. I. Menshikov and L. I. Ponomarev, Phys. Lett. B 167, 141 (1986). 

[11] L. I. Menshikov, Fiz. Elcm. Chastits At. Yadra 19, 1349 (1988), [Sov. J. Part. Nucl. 19, 583 
(1988)]. 

[12] J. S. Cohen and M. Leon, Phys. Rev. A 39, 946 (1989). 

[13] M. P. Faifman, L. I. Menshikov, and T. A. Strizh, Muon Catal. Fusion 4, 1 (1989). 
[14] M. P. Faifman and L. I. Ponomarev, Phys. Lett. B 265, 201 (1991). 
[15] M. Leon, Phys. Rev. A 49, 4438 (1994). 

[16] Yu. V. Petrov and V. Yu. Petrov, Phys. Lett. B 378, 1 (1996). 

[17] Y. Averin, V. R. Bom, A. M. Demin, et al., Hyp. Interact. 138, 249 (2001). 



30 



[18] N. Kawamura, K. Nagamine, T. Matsuzaki, et al., Phys. Rev. Lett. 90, 043401 (2003). 

[19] M. C. Fujiwara, Ph.D. thesis, University of British Columbia (1999). 

[20] M. C. Fujiwara, A. Adamczak, J. M. Bailey, et al., Phys. Rev. Lett. 85, 1642 (2000). 

[21] T. A. PorcelH, A. Adamczak, J. M. Bailey, et al., Phys. Rev. Lett. 86, 3763 (2001). 

[22] G. M. Marshall, A. Adamczak, J. M. Bailey, et al.. Hyp. Interact. 138, 203 (2001). 

[23] W. E. Lamb, Phys. Rev. 55, 190 (1939). 

[24] K. S. Singwi and A. Sjolander, Phys. Rev. 120, 1093 (1960). 

[25] L. Van Hove, Phys. Rev. 95, 249 (1954). 

[26] K. Fukushima, Phys. Rev. A 48, 4130 (1993). 

[27] L F. Silvera, Rev. Mod. Phys. 52, 393 (1980), and references therein. 

[28] P. C. Souers, Hydrogen Properties for Fusion Energy (University of California Press, Berkeley, 
1986). 

[29] A. Adamczak and M. P. Faifman, Phys. Rev. A 64, 052705 (2001). 

[30] P. E. Knowles, G. A. Beer, G. R. Mason, et al., Phys. Rev. A 56, 1970 (1997). 

[31] L. L Menshikov, Muon Catal. Fusion 2, 273 (1988). 

[32] V. N. Ostrovski and V. L Ustimov, Zh. Eksp. Teor. Fiz. 79, 1228 (1980), [Sov. Phys. JETP 

52, 620 (1980)]. 

[33] N. T. Padial, J. S. Cohen, and R. B. Walker, Phys. Rev. A 37, 329 (1988). 
[34] U. Fano, Phys. Rev. 124, 1866 (1961). 

[35] C. Petitjean, D. V. Balin, W. H. Breunlich, et al.. Hyp. Interact. 118, 127 (1999). 
[36] C. Petitjean, Hyp. Interact. 138, 191 (2001). 

[37] S. E. Jones, S. E. Anderson, A. J. CafFrey, et al., Phys. Rev. Lett. 51, 1757 (1983). 
[38] M. Jeitler, W. H. Breunlich, M. Cargnelh, et al., Phys. Rev. A 51, 2881 (1995). 
[39] Yu. V. Petrov and V. Yu. Petrov, Zh. Eksp. Teor. Fiz. 100, 56 (1991), [Sov. Phys. JETP 73, 
29 (1991)]. 

[40] Yu. V. Petrov, V. Yu. Petrov, and H. H. Schmidt, Phys. Lett. B 331, 266 (1994). 

[41] A. Akhiezer and I. Pomeranchuk, Zh. Eksp. Teor. Fiz. 17, 769 (1947), [Sov. Phys. JETP 11, 

167 (1947)]. 
[42] G. C. Wick, Phys. Rev. 94, 1228 (1954). 

[43] S. W. Lovesey, Theory of Neutron Scattering from Condensed Matter (Clarendon Press, Ox- 
ford, 1984). 

[44] B. D. Josephson, Phys. Rev. Lett. 4, 341 (1960). 

[45] F. J. Mompean, M. Garcia-Hernandez, F. J. Bermejo, and S. M. Bennington, Phys. Rev. B 

54, 970 (1996). 
[46] H. Bethe and G. Placzek, Phys. Rev. 51, 450 (1937). 

[47] M. P. Faifman, T. A. Strizh, E. A. G. Armour, and M. R. Harston, Hyp. Interact. 101/102, 

179 (1996). 

[48] A. Toyoda, K. Ishida, K. Shimomura, et al., Phys. Rev. Lett. 90, 243401 (2003). 
[49] T. M. Huber, A. Adamczak, J. M. Bailey, et al.. Hyp. Interact. 118, 159 (1999). 
[50] A. Adamczak, M. P. Faifman, L. I. Ponomarev, et al.. Atomic Data and Nuclear Data Tables 

62, 255 (1996). 
[51] J. S. Cohen, Phys. Rev. A 34, 2719 (1986). 
[52] A. Adamczak, Hyp. Interact. 119, 23 (1999). 

[53] J. Wozniak, A. Adamczak, G. A. Beer, et al., Phys. Rev. A 68, 062502 (2003). 



31 



